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The use of techniques from loop quantum gravity for cosmological models may solve some difficult 
problems in quantum cosmology. The solutions under a number of circumstances have been well 
studied. We will analyse the behaviour of solutions in the closed model, focusing on the behaviour of 
a universe containing a massless scalar field. The asymptotic behaviour of the solutions is examined, 
and is used to determine requirements of the initial conditions. 
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■ I. INTRODUCTION 

o ■ 
o ; 

I On large scales, the description of the universe provided by general relativity is consistent with all measurements 
1^ . to date. The cosmological models that have been generated seem to account for much of the history of the universe. 
Q ' However, the cosmological singularity that appears in this approach has suggested that a quantum theory of cosmology 
is required. Most studies of quantum cosmology have involved mini-superspace models. Homogeneity (and sometimes 
isotropy) is used to limit the degrees of freedom of the classical framework to a finite number. The reduced phase 
0^ I space is then quantized. Unfortunately, most such quantizations have been unable to provide a consistent solution 
to the singularity and initial value problems. Furthermore, the absence of time has lead to additional interpretation 
problems. 

In hopes of resolving these issues, the tools of loop quantum gravity have been applied to cosmology in a similar way. 
Loop quantum cosmology (LQC) involves the same reduction of the phase space as other mini-superspace models; 
, however, the system is quantized using the same techniques as the full theory of loop quantum gravity. In other 
00 ' mini-superspace models, the quantization is not necessarily related to that of a quantum theory of gravity. This is 
', thought to be responsible for many of the difficulties that arise. 

■ The results claimed by loop quantum cosmology to date have been impressive. First of all, it has been argued 
' that the cosmological singularity is naturally removed in this framework [1]. Secondly, the initial conditions can be 

O imposed from the behaviour of the equation itself [2]. Finally, classically divergent quantities are bounded in this 
framework(e.g. [7]). This introduces modifications to the evolution of the universe that seem to produce inflation [3] 
and eliminate chaos in the Bianchi IX model [4]. 

However, much of the successes of LCQ have come from the semi-classical analysis of the results. The interpretation 
J> of the quantum states themselves is not entirely clear. In particular, as with the Wheeler-DeWitt equation, 'time' does 
" '~j not appear as an explicit variable in the quantum framework. In the flat model, one typically picks the scale factor to 
rN represent a 'time' variable. From this, one can consider the evolution of the other degrees of freedom with respect to 
that time. However, in closed models where the universe recollapses classically, this interpretation seems less intuitive. 
" " In particular, the same volume corresponds, in general, to two different times. Furthermore, the classical solution 
attains a maximum volume, which one would expect to apply in the quantum case as well. Therefore, the probability 
should decay for large volumes. 

Much of the work to date in LQC has involved fiat models. This is not unreasonable since measurements indicate 
that we will in a very nearly fiat universe. However, if one wants to address current cosmological models, one must 
be capable of addressing the 'fiatness problem' associated with big bang cosmologies. As such, one must be able to 
understand the fiat model as a limit of the closed model. 

Because of the many issues involving closed cosmological models, it seems necessary to understand the nature of 
these solutions in loop quantum cosmology. In this paper, we will examine the solutions of the closed model in loop 
quantum cosmology in the presence of a scalar field. In particular, we will focus on the asymptotic behaviour of the 
model and how one imposes initial conditions in the exact solutions. Our analysis suggests that one cannot impose 
initial conditions that produce bounded solutions without introducing additional problems. 



II. THE CLOSED ISOTROPIC MODEL 



Loop quantum gravity is constructed by quantizing the phase space generated by the three-space connection variable 
A and its conjugate E. In the isotropic model, these variables can be written in terms of the symmetry generators of 
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the space and real multipliers. The end result is a phase space defined by the real variables c and p [7] such that 



{c,p} = ^K-f. (1) 

where k = SttG and 7 is the Barbero-Immirzi parameter, c represents the degrees of freedom associated with the 
connection A and is given by, 

c = T-^K, (2) 

where F is the spin connection and K is the extrinsic curvature. The spin connection is related to the intrinsic 
curvature of the spacial hypersurface, such that F = 0, 5 in the flat and closed models respectively. The conjugate 
momenta, p, is related to the scale factor a through \p\ = a^. With the Hamiltonian constraint, 

H = -±-sgn{p)^\ic' - c + = (3) 

this classical phase space describes isotropic cosmologies in general relativity. 

When quantizing the phase space of the full theory, one uses smeared variables rather than the basic variables 
themselves [5]. Since one would like to follow the full theory as much as possible, the holonomies of the reduced phase 
space can be defined in a similar way. In this case, the integration can be taken over straight lines because of the 
symmetries, giving the holonomies of A [7] as 

he = PexpiJ^A) = cos(cm) + 2sin(c/i)(e"a;^)r' (4) 

where /x e (—00, 00) gives the length of the integration (which is not a physical length) and P denotes path ordering. 
This leads to an orthonormal basis of states defined by 

|M)=exp(^). (5) 

From here, one defines p = —i^j^ which is an eigenfunction of these states with eigenvalue = — g^- From the 
classical relation V = \p\^^'^, one can define 

V^lM) = (^)^/^|M)-V.lM). (6) 

This is a well-defined operator since p is diagonal. In fact, this can be exploited to produce an operator for any 

positive power of p. From these simple operators and the holonomies, one can produce operator representations 
of more complicated quantities. The action of holonomies on the basis states can be determined by expanding in 
imaginary exponentials and using the relation 

e*«|M) = |M+l). (7) 
Using the above relations, one can define an inverse volume operator 

fc — —j 

where j and I parameterize the ambiguities involved in defining the operator from classical expressions [11,12], 

where I € (0, 1) and j is a positive integer. The integer j arises from out freedom to use different representations of 

SU(2). The parameter I arises from the fact that classically = However, the operators for will 

not preserve this identity on all scales. Asymptotically, this expression will behave like V~^, for all allowed values of 
I and j, since 

p',+m-pi-n^ = i^Yia + - (1 - ^Y] - (^)'(2Mm'-^ (9) 
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when II m. However, on very small scales, the inverse volume operator is bounded where the classical value is not. 
The precise behaviour of the eigenvalues on these small scales depends critically on j and /. 

In order to determine the physical states and obscrvablcs we need to define a Hamiltonian constraint operator. One 
can try to construct the operator along the lines of the full theory [5] where the Euclidian and Lorentzian constraints 
are quantized separately. This procedure was performed in [9] for both the flat and closed models. For the flat model, 
this approach works well; however, in the closed model, the quantizations is unstable [10]. Therefore another approach 
was found to be necessary for the closed model. 

There is a different approach is based on the techniques used for the Bianchi IX model [14,15]. Since F = i for the 
closed model, one can produce an operator for the Hamiltonian constraint using c — F in place of c. This method is 
likely to have less in common with the full theory since it exploits the form of F, but we can use the same techniques 
as the flat case. The resulting Hamiltonian constraint operator was found to be [8] 

C,rar = 96i{j' t^Hl)-' {sbi' {t^o^^^) cos' {f^o^-^^) ' ^(F^ - F))(sin(Mo|) V cos(mo|) - cos(Mo|)Vsin(Mo|)) (10) 

where the sines and cosines have the same action on states as previously described, fig G (0, oo) is an ambiguity 
parameter related to the 'length' of holonomies used to construct the constraint. The F terms in the trigonometric 
operators will be pulled out to form a complex exponential multiplier in the eigenvalue. As in the full theory, any 
matter Hamiltonian can be added directed to this operator to form a full Hamiltonian constraint. 

From the constraint operator proposed in (10), one can define the physical states of the system; ie, those that are 
annihilated by the constraint operator. The obvious starting point is the volume eigenfunctions. The action of the 
constraint operator on a volume state is 

Cgravln) = 3{'y\ll)-\V^+^^ - V^.^Jie''^^^ \n + A^o) - (2 - ^I^hHt' - F))|m) + €-''^^"1^ - 4^0))- (11) 

Clearly the volume eigenstates, other than the zero volume state, are not annihilated. Therefore, we will write any 
physical state as a sum over volume states via 

l*) = E^''WI'")' (12) 
ft 

where </> denotes the matter dependence. Therefore, imposing C\^) = yields the condition 

- WJe-'2r'*<'s^+4^„ - (2 - 4M^7'(r' - r))(W - Vi^-i^o>^ + ( Vs^o - V^.,^Je'^^^^ s^.^^^ 

= —l^h^KllHmauS^,, (13) 

where the matter Hamiltonian, H^att^ has been added directly to the constraint, as in the full theory [6]. This is 
the discrete equation that describes the state of the universe. Therefore, given a matter Hamiltonian, we use (13) 
to determine the coefficients of the physical state. However, one should keep in mind that the universe is not in one 
volume state that 'evolves' according to this equation. The states are only those that include the sum of all volume 
states. In the flat model, the parameter fj, is frequently used to represent 'time'. However, it is not clear that this 
interpretation applies in general. 



III. EXACT SOLUTIONS 



In order to find exact solutions, we will need to choose a form of matter. For simplicity, we will use a scalar field 
with the Hamiltonian 

Hm.u = \a-hl = ~a-^^„ (14) 

where o~'^ is defined as in (8). We assume that the coefficients have the form 

s^{4>) = c^e'-^ (15) 

where lu is a constant and c is independent of (f). This particular form of the matter Hamiltonian and the intial wave 
function were chosen to make it possible to find a solution. Nevertheless, we can write any possible matter dependence 
for this Hamiltonian as sum of solutions of the form (15). Using this form, we can rewrite the RHS of (13) as 
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- ^^J'h^KllHmattS^, = -lll^Hlu^dMSi,, (16) 

where 

.^ = 5^. (IV, 

We also introduced an exponential containing T to eliminate the exponentials from the Hamiltonian constraint 
equation via the transformation s — > se^'^s . The resulting equation can be written in iterative form 



5^+4^„ = [(2 - ^H„l{T^ - r))( W W )^^ _ ^^^isgn{i,)>^HldM^, 



-( 73 V^^j 



which holds for all ji ^ —^z^io- None of these equations depend on sq which drops out of the equation entirely. This 
is why we have no equation for ji = —Afic,. Nevertheless, one can evolve though the classical singularity, ^ — Q, since 
choosing /it = in (18) is simply 54^^ = s_4^^. In solutions where /i ^ 4/Lto, (18) is always valid; thus, these are well 
defined around the zero volume state. The sgn{p) term is introduced because one has the freedom to choose the sign 
on the Hamiltonian constraint. The sgn{ix) was chosen since there is no reason to expect different solutions in the 
classically disjoint regions a < and a > 0. 




FIG. 1. This is the exact solution to the closed F = | with w = 55 and j = 20. The initial value was set a = Afio with 
Ho = Vs/A. The state clearly diverges as we go to large scales. Since large p is classically disallowed, we would like not to 
admit these states as physical possibilities. The factor exp(i/ur/2) has been left out of the solution. 
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One can solve this equation itcratively by choosing 'initial' conditions at any particular point. In this model, we 
have 2 degrees of freedom, as long as the solution satisfies the equation at /x = 4/Zo. However, the form of the solutions 
shows the first problem with this approach. We have factored out the e'''^ which will be present in all solutions. 
This factor implies that the state changes rapidly on scales near the plank scale. These oscillations on small scales 
also appear in the flat model and remain to be understood. However, in the flat model, one avoids such solutions 
on the grounds that they do not produce semi-classical physics. In this case, since the terms appear explicitly in the 
equation, we can still approximate (13) with a differential equation on large scales. Therefore, this is only a problem 
in the sense that it is not clear what this phase represents and how it relates to quantities that can be 'measured'. 

Before we get to the properties of the exact solutions, we will look at the approximate behaviour on large scales. 
This will give a much better qualitative understanding of the solution than can be seen from the difference equation 
itself. 

First of all, let us define the new variable such that 



t^. = e"'^V^+^^-V^.^Js^. (19) 



This satisfies the simpler equation 



W4^o - + V-4^o = -^^t^llldM— - ^H{T^ - r)fj,X. (20) 

I V+z^o ''n-tJ.o ) 

When /X is large, the LHS of (20) can be approximated by a differential operator, 



- + *^-4m„ = 16Mo0 = ^(^^)''"o0- (21) 



In the same limit, we have dj^i^jj) = p The final approximation for the difference in the volume eigenvalues. We 
can rewrite these in terms of p using, 



W - = (^)^/^((i + ^r^' - (1 - « C^f/'s^t^o (22) 



and y/JI = ,/-^. Prom here we get 



W-W = 4^- (23) 



The end result is the Wheeler-DeWitt-like equation. 

In the fiat case (F = 0), one can find a solution by assuming t = Ap^. The equation gives the solution 



From this solution, one sees that there is a critical value of w, ui'^ = above which we get oscillatory solutions. 
Assuming co > ujc, the solution takes the form 

t = y/p{B exp{+iRlog{p)) + Cexp(-ii?log(p))), (26) 

where R is the imaginary part of r. Recall that we absorbed volume terms into t. From (23), we find that s cx 
Thus, in the large-scale limit, s is purely oscillatory. Similar analysis has been performed for the fiat model using 
quantizations of the Euclidian [16] and the Lorentzian [13] constraints and yield the same critical value. 

In the closed case (F = i), one can drop the matter term for large p since it drops off as relative to the other 
terms. The resulting equation has simple solutions in terms of exponentials, 

tin) = Bexp(- Ji^ct) + Cexp(+j7/x) (27) 
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where B and C are constants. The real exponentials are not surprising, given that large values of a are classically 
disallowed. Like any conventional quantum mechanics problem, one would like to impose the boundary condition that 
f — > as 1^1 ^ oo. Clearly one would set C = for large positive fi, and S = for negative /U. However, the iterative 
equation defines solutions for both positive and negative ji. 

Let us consider a solution where we pick some initial value 54^^ = De^'^'^. From (18), we find that 



s,,„ = [(2 + /'""n - l'lAs9n{ij)u^HldMW^o- (28) 



Since Sg/i,, (x D, D will scale all values of the solution for positive fi. Furthermore, when we set /U = 0, the constraint 
equation turns into 54^^ = The evolution is symmetric in /z so f will be present for all negative /x as well. 

Therefore, the solution will not depend on how we choose our initial value. One such exact solution is shown in figure 
1. As one can see, the solution diverges at large fi. Since we can only set an overall normalization, the only solution 
we have (other than \^) =0) is unphysical. 

In the case where we define our initial value at fi ^ nAiio, where n is an integer, there is no additional constraint on 
the solution from the zero volume state. However, because we no longer have the zero volume to eliminate terms, we 
do not have symmetric evolution in positive and negative /i. Therefore, if we eliminate the clivcrgcncc in the negative 
domain, we will still have the problem in the positive domain and vice versa. Figure 2 shows the state of the universe 
when we impose that — > for negative fj,. Clearly, the problem remains. This is very similar to the problem of 
producing the pre-classical solution in the flat model with the Lorentzian constraint [13]. In that case, one can only 
produce smooth solutions in either /z < or /U > 0, but not both simultaneously. 
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FIG. 2. This is the exact solution to the closed model F = 1/2 with a; = 55 and j = 20. The initial value was set a = 3/Uo 
with /Xo = 1. The initial conditions have chosen to eliminate the growing solution in the —fj, domain. However, it is clear we 

cannot eliminate the growing solutions in both positive and negative domains for the asymmetric model (i.e. the model where 
we don't hit zero volume). The factor exp(i/ir/2) has been left out of the solution. 
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Solutions to this particular model have been presented [8] for the case where a cosmological constant is included. 
The solution avoids the difficulties presented here. In particular, one can rewrite the matter Hamiltonian as 



(29) 



where Vo is the cosmological constant. The divergent solutions described above are not observed in this case because 
exponentially growing solutions are no longer a feature of the asymptotic solution. In particular, this introduces an 
additional term which dominates at large ii. The resulting behaviour is described by 







where V 



Vo-^l^ol^i^ip- The solutions to (30) are defined in terms of Airy functions 

A'Ai{-V^/^\fi\) + B'Bi{-V'/%i\), 



t : 



(30) 



(31) 



where A',B' are constants. For Vo > 0, the solutions are well behaved for all /x. Coincidently (or perhaps not), 
this is the same condition that determines if large a is classically disallowed. For a generic matter Hamiltonian and 
wave-function (j) and a will not decouple, but there will be regions of the a — (p plane that are classically disallowed, 
where we expect similar problems will arise. 

There is reason to be concerned that the numerical results presented in figures 1 and 2 do not represent the exact 
solutions of the constraint equation. Because the solutions will take the form of exponentially growing and decaying 
solutions, numerical round-off errors may grow exponentially. This problem was tested explicitly. First of all, the 
solutions calculated using between 16 and 1000 digits all behaved in a similar manner and were within a few percent of 
each other where the divergences become large. Secondly, one can produce solutions that are decaying where solutions 
usually diverge by avoiding the constraint at zero (or leaving the opposite /x domain to diverge). These solutions are 
qualitatively different and near fj. = the ratio s^^^/s^^^ differ by orders of magnitude. This behaviour has been 
tested over a reasonable range in 7, jig, I, j and u with similar behaviour in all cases. 
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FIG. 3. This is the exact solution to the closed model F = 1/2 with uj = 55, j = 20 and 7 = ^^l;^- The crosses represent the 
solution for I = .45 and solid line is for / = .4475. The initial value is set sX n = Ajio and is equal in both cases. The change in 
the sign of the growing exponential part of the solution suggests that we can fix / such that only the decaying solution remains. 



Although these ambiguity parameters provide additional freedom in our solutions, it is unclear whether they can 
cure our problems or whether they simply relocate them. First of all, there have been arguments presented that should 
fix our choice of both 7 and /Xq [7]. Assuming these to be fixed, the only remaining parameters are j and I. If we fix 
all values and adjust I over the allowed range (0, 1), there are changes in the sign of the diverging exponential for fixed 
S4^^ . This suggests that we can fine tune I such that the exponentially growing solution vanishes. Unfortunately, this 
fine-tuned value depends on lo. 

In figure 3, we have an example of the sign flip that occurs for a very small change in I when w = 55. However, 
if a; = 50, the solutions for the same two values of I are virtually indistinguishable, as seen in figure 4. Further 
choices of I seem to suggest that there is no common I in that range for both = 50 and = 55 that eliminates the 
unbounded solution. The same holds for a wide range of choices of w < 500, j < 1000 and 7 < 1. Since the inverse 
density operator is defined on the kinematic Hilbert space, its definition should not depend on the form of the matter 
Hamiltonian or the particular state of the matter. In this particular model, one should be able to expand the solution 
as a sum over all lo. One should not have to define l{w), before doing so. Therefore, it seems unlikely that a choice 
of these parameters produce an acceptable model. 

The bounded solutions achieved by fine-tuning I have the feature of being sharply peaked around /x = 0, as suggested 
by figure 5. This is an additional difficulty of the fine-tuned solutions. Not only does one need to fix 1{lu), but one also 
needs to explain to meaning of the solutions. In particular, one should address why the only states are dominated by 
the smallest volume eigenstates. 
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FIG. 4. This is the exact solution to the closed model F = 1/2 with uj = 50, j = 20 and 7 = ^^|^- The crosses represent the 
solution for I = .45 and solid line is for / = .4475. The initial value is set at = 4/Xo and is equal in both cases. These solutions 
are nearly identical for all //. This suggests that there is no single value of I that will produce bounded solutions for all choices 

of CJ. 

Beyond fixing I, we could also fix other parameters in hopes of producing bounded solutions. For example, one 

could adjust the initial point (instead of assuming the symmetric solution) to some critical fi G xAfio where x G (0, 1). 
This would eliminate the continuous spectrum of solutions for fixed uj, etc in the same way we fixed I or any of the 
ambiguity parameters. However, any one of these choices would depend on all the other parameters in the theory. 
Therefore, we can restrict one variable in order to produce bounded solutions but at the cost of introducing additional 
matter and ambiguity dependences into the states and the constraint. Therefore, one needs to ask whether it is 
sensible to have such a theory. 



o 
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FIG. 5. This is a plot of the log of modulus of for the closed model F = 1/2 with cj = 55, j = 20 and 7 = The 
crosses represent the critical solution for I = .447.... in which the solution is bounded at large scales. The state is normalized 
such that the amplitude of the oscilation is of the same order as the other solutions. Thus, the value at yu = 4/io is nearly 10®° 
compared to 10~®° in the other solutions shown. Clearly, the behaviour near fj. = dominates the solution. 



IV. DISCUSSION AND CONCLUSION 

In the analysis of the behaviour of a massless scalar field in the closed model from ILQC, it seems that divergent 
behaviour of the solutions at large scales is unavoidable. In this analysis, we consider the solutions for a wide range 
of parameters (e.g. j, I). In every case, the solution oscillates at small scales but grows exponentially for large /i. 
Attempts to constrain the ambiguity parameters through the boundary conditions seems to require significant matter 
dependence, and thus is an undesirable solution, if it works at all. 
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When trying to produce a closed model using techniques from the full theory, one gets an unstable discrete equation. 
If the symmetry is used to make the constraint simpler to quantize, we get unbounded coefficients for the volume states 
in classically disallowed regions. If we are to interpret these solutions as representing a probability, then clearly we 
cannot allow the latter quantization. This assumption is not completely trivial. The physical meaning of these states 
is still unknown. Without an inner product and observables, the interpretation of the solution is not entirely clear 
(one possible interpretation can be found in [17]). Nevertheless, it seems unlikely that one can find an interpretation 
where the large-scale divergences are appropriate. 

We are left with the problem of how to quantize the closed model in loop quantum gravity. From this point, it 
seems unlikely that a minor modification to the techniques applied to the fiat model will be successful. Therefore, 
one needs to reconsider the techniques used from the first step in this procedure. One of the problems is that the 
mini-superspace treatment has removed the coordinate freedom one usually uses to eliminate the spin connection. If 
this is indeed the case, there must be a way to construct the cosmology from the full theory that does not suffer from 
this problem. Hence, the full theory might avoid these difficulties. This would also suggest that this mini-superspace 
model has less in common with the full theory than was intended. 

Perhaps the most interesting aspect of the divergences that occur in this model is that they seem to occur under 
the same conditions that lead to a re-collapsing universe. The re-collapsing universe has typically been difficult to 
interpret in quantum cosmology. In this framework, it would be very difficult to understand the 'evolution' of a 
re-collapsing universe in terms of the volume states. It is interesting that our quantization seems to fail in the same 
place our framework would lead to significant interpretation problems. 

It remains to be seen whether these types of problems rise in the Bianchi IX model from which this isotropic 
quantization was derived [15]. The homogenous solutions contain more degrees of freedom but will also have additional 
constraints and boundary conditions. One might hope that the problems found here result from limiting the phase 
space to classically isotropic solutions. If the homogenous solutions avoid these divergences then it would suggest that 
we might not need to worry about the difficulties in the isotropic case. 
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